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Abstract 

Consider the Jacobi operators J given by {Jy)n = anyn+i + bnyn + an-iUn-i, Un S 
(here yo = yp+i = 0), where 6„ = 6* and a„ : deta„ ^ are the sequences of 
m X m matrices, n = l,..,p. We study two cases: (i) a„ = a* > 0; (ii) o„ is a 
lower triangular matrix with real positive entries on the diagonal (the matrix J is 
(2?7i + l)-band mp x mp matrix with positive entries on the first and the last diagonals). 
The spectrum of J" is a finite sequence of real eigenvalues Ai < ... < Xjy, where 
each eigenvalue Xj has multiplicity kj ^ m. We show that the mapping (a, b) i— > 
{(Aj, kj)}i © {additional spectral data} is 1-to-l and onto. In both cases (i), (ii), we 
give the complete solution of the inverse problem. 



1 Introduction and main results 



We consider the finite matrix- valued Jacobi operator J' = J'a,b acting in (C^)™ and given by 



( h\ a\ 
a\ h-i a2 
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bp-1 



where b = 

bneS = {b 

(i) a„ G S+ 

(ii) a„, G Li 



. 

y ... w w u,p_^ 

an)i~^ are the finite sequences of m x m matrices such that 
b = b*} and det a„ 7^ for all n e l,p = {1, 2, . . . ,p}. We consider two cases: 
- {a : a = a* > 0}; 

= {a : a is a lower triangular matrix with real positive entries on the diagonal}. 
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In both cases we obtain the complete characterization of the set of spectral data that 
correspond to these classes of operators Jaf,- Note that in the second case Ja,h are self-adjoint 
(2m-|-l)-band matrices of the size mp x mp with positive entries on the first and the last 
diagonals. 

The spectrum a{J) oi J — J* \s, finite sequence of real eigenvalues 

Ai < A2 < . . . < Ajv, 

where each eigenvalue Aj, j e 1, A/" has multiplicity kj e l,m, i.e., kj is the number of the 
linearly independent eigenvectors corresponding to A^. Note that 

kx ^ k-i ^ . . . ^ kf^ = mp. 

Introduce the fundamental m x m matrix-valued solutions ip{z) = {fn{z))V~^ said xi^) = 
{Xn{z))Q^^ of the equation 

anVn+l + KVn + = ^^n, n) G C X T~p, (1.2) 

such that 

<^o = Xp+1 = 0, ^i = Xp = I, 

where / is the identity m x m matrix and we set = ap = I for convenience. Note 
that i^niz) and Xn('2), ri G l,p are matrix-valued polynomials such that degLpn{z) = n, 
degXn(-2) = p+l — n. Each eigenvector ■0 = {ipn)i, J'4' = ^ji^^ has the form ipn = '~Pn{^j)v 
for some v e C"*. Hence, the eigenvalues of J are zeros of det 

Definition 1.1 (Spectral data). For each eigenvalue Xj, j G 1, A^, we define the subspace 

Sj = Ker <^p+i{Xj) = |/i G C"' : ipp+i{Xj)h = o| C C"", dimSj = kj ^ m, (1.3) 

the orthogonal projector Pj : — > Ej C onto £j and the positive self-adjoint 
operator gj : £j — > £j given by 

p 

9j^Gj\£., where Gj ^ Pj^^l{Xj)^n{\)Py (1.4) 

n=l 

We now describe the connection between our spectral data and the matrix-valued Weyl- 
Titchmarsh function M[z) given by 

M{z) = -Xi{z){xo{z))-\ zee. (1.5) 

Proposition 1.2. The function M{z) — M*{z) is analytic inC \ {Aj,j G l,N}. Moreover, 

TV N 

^(^)--E7zV' E^^ = ^' (1-6) 

1 Z Aj . 
j=l J j=l 

where the self-adjoint matrices Bj = res M{z) = B* are given by 

z=Xj 

Bj\e,=9i\ i?.lc-e£, =0, J e MV. (1.7) 
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In order to formulate our main result we need 



Definition 1.3. We call the system {{Xj,Pj),j G 1, A^} of the distinct real numbers Xj and 
the orthogonal projectors Pj : C"* the p — tame system, z/^^^^rankPj = mp and 



det 



/ To Ti ... T,., \ 
Ti T2 ... Tp 



N 



^0, Ts = J2^jPj^ s = 0,..,2p-2. 



\ 7p-i Tp ... T2P-2 I 

Remark. (i) The mp x mp matrix (Tsfc)^'^^Q is always non-negative definite (see p.ip ). 
Hence, the system {(Aj, Pj), j G 1, A^} is p-tame iff this matrix is strictly positive definite. 

(ii) Let N = p and kj = dim£^j = rankPj = m for all j G l,p, i.e., Pj = I for all j. Then for 
each distinct values Ai < . . . < Ap the system {{Xj, G l,p} is p-tame. 

(iii) Let p = 2, = 3, ^1 + ^2 = ^3 = m. Then the system {(Ai, Pi), (A2, P2), (A3, /)} is 
2 -tame iff Ker Pi n Ker P2 = {0}. 

Introduce the set of spectral data 

{p ^ N ^ pm, Ai < . . . < Aat are real numbers; 
p xjv . {(Aj,P,), j G 1,A^} is a p- tame system; 
V i' j^yjJi ■ g. ■ RanPj — > RanPj- are linear operators such 
that gj = g*>0 and Zf=i P,9j'P, = I 

and the mapping 

^:((a„)r;(M?)--(A„P„^7,)f • 
We formulate our main result. 

Theorem 1.4. (i) The mapping ^ : S^"' xSP ^Sp is one-to-one and onto, 
(ii) The mapping \I' : L^^^ x ^ Sp is one-to-one and onto. 

There is an enormous literature on inverse spectral problems for scalar (i.e., m = 1) 
Jacobi matrices (see book [Tj and references therein), but considerably less for matrix- 
valued Jacobi operators (see ICXjRj and references therein). The inverse problems for finite 
scalar Jacobi matrices were considered by several authors (see |GSj and references therein). 
Some uniqueness results for matrix- valued Jacobi operators were obtained in |CGRj . and 
the intimate connection to matrix-valued orthogonal polynomials and the moment problem 
were treated in [C]. |DLj . 

The main goal of our paper is to give the complete characterization of the set of spectral 
data for operators J7a,6. We hope to use similar ideas to obtain the complete characterization 
of the spectral data for the Sturm-Liouville operators TCy = —y" + V{x)y, y(0) = y{l) = 0, 
on the unit interval [0, 1], where V = V* is a. m x m matrix potential. Note that the "local 
characterization" of the spectral data for operators Ti was obtained in |("Kj . 

In the proof we use the approach from [(jSj and some technique from |(yKj . The in- 
verse spectral problem consists of the following parts: i) uniqueness, ii) characterization, 
iii) reconstruction. In Theorem 11.41 we solve all these problems simultaneously. 
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2 Preliminaries 

For each two sequences of polynomials 'd{z) = (??„(z))q^"^, r] = {Vniz))^^^ we introduce the 
Wronskian 

{^,v}n{z) = i^*^{z)anr]n+i{z) - i9*^^-^{z)alrin{z) , n = 0, ..,p. 
If ()1.2|) holds for both and t], then {■!?, r]}n{z) does not depend on n. In particular, 

^*p{z)^P+iiz) - '^l^-^{z)^p{z) = {v?, i^}{z) = and Xo(^) = {x, = fp+i{z), (2.1) 

since oq = Op = / by our convention. Recall that 

Sj = Ker ^p+i(A,) = {/i G : ^p+i{\j)h = 0}, j G i;iV, 
and Pj : — is the orthogonal projector. Also, we need the subspaces 
= Ker(^;+i(A,) = {/i G C'" : y.;+i(A,)/i = 0}, j G T;1V, 

and the orthogonal projectors Pj : £j. Note that dimf"] = dim£^j = kj. 

Lemma 2.1. For each j &1,N the following identities are fulfilled: 

Py{X,)Pj = Vp{\)Pj. Xi{\)Vp{\)Pi = Pj, (2.2) 
p 

^ ^Yl '^ni^j)'^n{Xj) = <^*p{Xj)iPp+l{Xj) - '^*p+l{Xj)^p{Xj). (2.3) 

Proof. Using (|2.1|) . we obtain ^Pp_^_^{Xj)^Pp{Xj)Pj = ip*{Xj)ipp+i{Xj)Pj = 0. This means 
^p{Xj)Pj C S^j, i.e., P^ipp{Xj)Pj = ipp{Xj)Pj. Let 

Vn = XniXj)(Pp{Xj)Pj - (pniXj)Pj, U = 0, ..,p+l. 

The sequence {yn)o^^ satisfies fll.2|) for z = Xj and 

Vp+i = -(Pp+iiXj)Pj = 0, yp = ipp{Xj)Pj - ipp{Xj)Pj = 0. 

This yields y„ = for all n G In particular, yi = xi{Xj)(Pp{Xj)Pj — Pj = 0. 
Furthermore, using equation ()1.2p . we obtain 

f*n+l{z)an + V^n(^)(^n " ^) + <^;_i(^)a„_i = 0, 
an^Pn+l{z) + {hn - z)ipn{z) + a*^_^<pn-l{z) = (Pn{z) . 

Multiplying the first equation by (pn{z) from the right, the second equation by >fn{z) from 
the left and taking the difference, we deduce that 

^*niz)^niz) = {(p,ip}niz) - {(f , (p} n-liz) , U G l,p. 

The summing implies fl2.3p since {ip, 0}o = 0. Note that / = ipl{Xj)ipi{Xj) by definition. □ 
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Lemma 2.2. (i) For each j E I, N the mapping Yj = Pj ipp^i{Xj) Pj : Sj £j is invertible. 
(a) For each j E 1,N the following asymptotics holds: 

PYr^P^ 

(V9p+i(2))-^ = ' \ ' +0(1) as z^\. (2.4) 
Z - Aj 

Proof, (i) Let Yjh = for some h ^ Sj. Using ()2.3|) and ()2.2|) . we obtain 

p 

n=l 

Since the left hand side is positive definite on Sj, the operator Yj = Pjipp+i{\j)Pj is invertible. 
(ii) Let 

" A{z) B{z) 
C{z) D{z) 



where 



C = P^^p+.Pj- : Sf^S], D = P]^p+,P,: S,^S]. 



Since ipp+i{\j)Pj = and Pj(^p+i{\j) = 0, we have B,C = 0{t) as t = z — Xj — > 0. Moreover, 

D = t- Pl^p+i{\j)P, + 0{e) = tY, + O(t'), A = X, + 0{t), 

where Xj = {Pj)-^ipp^i{\j)Pj- = ipp^i{Xj)Pj- : Sj- — > {^j)'^- The operator Xj is invertible, 
since i^p+ih = implies h G The operator Yj is invertible due to (i). Therefore, the 
standard formula for the inverse matrix 





B ^ 






D J 





A-^ + A-^BH-^CA-^ -A~^BH-^ 
-H-^CA-^ H-^ 



H = D- CA-^B, 



gives 

i^M^)) = [ 'o{l) t-%r^ + 0{l) ) t-.0. 
In particular, res {(fp+i{z))~^ = PjYf-^Pj. □ 

z=Xj 

Proof of Proposition 11.21 Note that Xo(j)xi{z) - Xii'^)Xo{z) = {x,x}o{z) = 0. This 
gives M*{z) = M{z), z e C. Due to dSl), the function (xo(^))^^ = {^p+i{z))~'^ has a 
simple pole at each point z = Xj. Therefore, the function M{z) = —Xi{z){Xo{z))~^ has a 
simple pole at z = Xj and 



B, = - r_es M{z) = xi(A,) ■ ( r_es (ifp+iiz))-')* = Xi{Xj) P^ P, 

Z At z — A ■)■ 
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where Yj = Pj (pp^i{Xj) Pj : £j £j. This yields BjI^nq^. = 0. Furthermore, using ()2.3|) and 
()2.2|) . we obtain 

G, = G* = P,^;^,{X,)^p{X,)P, = P,0;^i(A,)P,Vp(A,)P, = PjY*P^ . MX,)P,. 
Hence ((2I2D yields, 

BjG, = xi(A,)P,Vp(A,)P, = Xi(A,)^p(A,)P, = P„ 

i.e., -Bjl^;. = fi'j""'^. Asymptotics M{z) = —Iz^^ + 0{z^^) as z ^ oo (see (II. 2j) for n = 0) and 
the standard arguments from the function theory give ()l.fi|l . □ 



3 Proof of Theorem 11.41 



Lemma 3.1. The system {(Aj, Pj),j G 1, A^} isp - tame, iff there is no nonzero vector-valued 
polynomial F{z) G C" such that deg F ^ p—1 and PjF{Xj) = for all j G 1, N. 



Proof. Note that for each vector v 



)n ^, t's G C™", we have 



(^0 



p-1 



/ ^0 



p-1 



Ti 



T 



^2p-2 / 



^^0 
^^1 



V / 



N 



N 



= E < E ^^^'^-'■^'^ = E ^*(A,)P,F(A,) ^ 

s,k=0 j=l j=l 

where F{z) 

F{z) of degree at most p—1 such that PjF{Xj 



(3.1) 



Yll=o^''^^s- Therefore, det(rs+fc)s fc=o 7^ iff there is no nonzero polynomial 



for all j G 1, A. 



□ 



Proof of Theorem 11.41 The proof is similar for both cases (i),(ii) and consists of three 
parts. We need also two simple technical Lemmas 13. 2| 13.31 that are located at the end. 

1. * maps S^^x SP and L^"^ x Sp into Sp. 

Recall that the identity PjQj^Pj = J2f=i Bj = is proved in Proposition 11.21 In order 

to check that {(Aj, Pj),j G 1, A^} is p-tame, we use Lemma f3. II Suppose that PjF{Xj) = 
for all j G 1, A and some vector- valued polynomial F{z) of degree at most p — 1. Using 
asymptotics ()2.4p near poles of {xo{z))~^ = {'^p+i{z))~^ , we deduce that the vector- valued 



function (xo(^)) ^F{z) is entire. On the other hand, Xo'^{z)F{z) = 0{z~ 



p 



z- 



p-1 



as z oo. The Liouville Theorem gives Xo '^i^)F{z) = for all z G C and hence F 
2. * : X SP ^ 5n is one-to-one and * : LP^ x ^ S„ is one-to-one. 



) = 0{z-^) 
0. 



Following jGSj . we introduce the sequence of M-functions 



-Xn(2)[<_iXn-l(2)] 



-1 



n G l,p+l. 
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Using ()1.2p . it is easy to see that 

- (M„(2;))"^ = [{z - hn)Xn{z) - anXn+i{z)]{Xn{z)Y^ = I z - + a„M„+i(2;X (3.2) 

for all n G l,p. In particular, Mn{z) = —z~^ + 0{z~'^) as, z ^ oo for all n. Due to 
Proposition II. 2[ the spectral data {Xj, Pj, gj)i uniquely determine the function M = Mq. 
Therefore, the matrices bi and A = aia^ are uniquely determined by the asymptotics 

M{z) = -Iz~^ - biz~^ - {aia[ + &i)^~^ + 0{z~^) as z ^ oc. 

If tti, ai G 8+ (the case (i)) and aia* = A = aial, then ai = oi. The same is true, if ai, oi G 
L+ (the case (ii), see Lemma ESI)- Hence, the matrices bi and ai are uniquely determined 
by the spectral data. Therefore, the function Mi(z) = ai^{Iz — bi — (M(z))~-^)(a^^)* is 
uniquely determined too. Repeating this procedure, one uniquely determines all matrices 
bn, ttn in both cases (i), (ii). 

3. * : X SP ^ 5p is onto and * : L^^ x ^ Sp is onto. 

Let Pj, gj)i G Sp. We shall construct 6„ and a„ step by step. Introduce the function 

N R 

M{z) = -J2^, (3.3) 

j=i ^ 

where Bj = PjQj^Pj = B* ^ 0. Note that M{z) = M*{z), z e C, and 

lmM{z) = —{M{z) - M*{z)) > for all Im^ > 0. 

In particular, det M{z) 7^ 0, if Im 2; > 0. This gives 

-lmM-\z) = {M-^y{z) ■ lmM{z) ■ M'^z) > 0, Im^ > 0. 

Moreover, -M-^{z) = Iz + 0(1) as 2; — * 00, since Ylj=i ~ ^- Therefore, the Herglotz 
representation theorem for rational matrix- valued functions (see Lemma f3.2|) yields 

^ n 

-M-\z) = Iz + C -^—^ (3.4) 

s=l ^ ~ 

-.TV 



for some C = C* and = ^ 0, s G 1, K. Denote J^, = e Ker and let 

Qs '■ — s> jFg c be the orthogonal projectors onto J^g. In accordance with ()3.2|) . we set 

bi = -C. 

Let Aj 7^ /is for all j, s. Then all poles of the meromorphic function det M{z) are {Aj}^ 
and all roots are {^is}i ■ Moreover, each A^ is a pole of the multiplicity dim£^j = rankPj and 



7 



each /is is a root of the multiphcity dimjF, = rankQs- Since det M{z) = {—z) ^ + 0{z ™ ^) 
as z oo, we have 

K N 

rank Qs = —m + rank Pj = m{p — 1). (3-5) 

s=l i=l 

If Aj = /is for some j, s, then the corresponding pole and the root (partially) compensate 
each other but (j3.5j) still holds true. 

If p = 1, then ()3.5|) yields K = and the reconstruction procedure stops. 

Let p > 1. We show that the system {{iJ's,Qs),s G l,K} is {p — l)-tame, this is the 
crucial point of the proof. Suppose that QsG{^s) = for all s and some vector-valued 
polynomial G{z) of degree at most p — 2. Due to Lemma EIH it is sufficient to prove that 
G = 0. In view of (|3.4p . the vector- valued function 

F{z) = M-\z)G{z) 

is entire and F{z) = 0{z-zP~'^) = O^z^^^) as 2; — > 00, i.e., F{z) is a vector-valued polynomial 
of degree at most p — 1. Since M{z)F{z) = G{z) is entire too, ()3.3|) yields PjF{Xj) = for 
all j G 1, A^. Due to Lemma f3. 11 this implies F = 0. Hence, G = 0. 

Define A = Ylf=i (note that ^4 = ^4* ^ 0). In order to show that ^ > 0, suppose that 
f*Af = for some constant vector / G C™, / 7^ 0. This gives Qsf = for all s E 1,K, 
which is a contradiction with Lemma lXTl Hence, A > 0. In accordance with (j3.2p . we choose 
oi such that 

aial = A. 

In both cases: (i) oi = \/A > and (ii) see Lemma ESI there exists a unique matrix solution 
Oi of this equation. Now we set 

M,{z) = -a-\h -Iz- {M{z))-'){a-'r = -f^ 7^ , 

Z Us 



where 



K 



Ds = a^^Ds{a^^y, sel,K, and ^ 



s=l 



Note that KerD^ = alKei Ds = a\J^^. Let Qs : ^ ) be the orthogonal 

projectors. Suppose that QsG{iis) = for all s E 1,K and some vector-valued polynomial 
G{z) of degree at most p-2. Then G{ns) G a^J^^, i.e., (5s[(a^)"^G(/Xs)] = for all s G 1,K. 
Due to Lemma this gives {al)~^G = and hence G = 0. Therefore, the new system 
{{fis,Qs),s G l,K} is (p — 1)- tame. Repeating the procedure given above, we reconstruct 
&2, 02, &3, 03 and so on. □ 
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Lemma 3.2. Let the rational matrix-valued function f{z) have only real poles and satisfies 
f(z) = f*{z), zeC, and lmf{z) = \. {f{z)-f*{z)) > 0, Imz > 0. Then f{z) has the form 

f{z) = D,z + C-Y,^^ 

for some K ^ 0, e C = C* and D^^ D*^0, s^O, .., K. 
Proof. Due to the identity f{z) — f*{z), we have 

j=l s=l j=l ^ 

for some D^ j = D*^ and C = C* . Using the condition Im/(2;) > as Im^; > near the 
points z = Xg and z = oo, we deduce that Dg j = 0, if j ^ 2, and Dg = —Dg i ^0. □ 

Lemma 3.3. Let A e S+. Then there exists unique matrix a e L+ such that aa* — A. 

Proof. Let 

, f B C* \ f b 

"^^[C D )^ ^-[c d 

where B,b are positive real numbers, C,c are vectors and D,d are (m — 1) x (m— 1) matrices. 
The equation aa* = A is equivalent to 

- - CC* 

bb = b^ = B, cb - C, dd* ^D-cc* ^D- 



B 

Note that D — CC* > 0, since A > 0. Now the problem is reduced to the similar problem 
for (m — 1) X (m — 1) matrices. The case m = 1 is trivial. □ 
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